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With SNO data [jl] on electron-neutrino flux from the sun, it is possible to 

OO 

derive the v e survival probability P ee (E) from existing experimental data of 

Super-Kamiokande, gallium experiments and Homestake. The combined data 

Q_i| of SNO and Super-Kamiokande provide boron v e flux and the total flux of all 

Qh' 

' active boron neutrinos, giving thus P ee {E) for boron neutrinos. The Homes- 

take detector, after subtraction of the signal from boron neutrinos, gives the 

•rH , 

flux of Be+CNO neutrinos, and P ee for the corresponding energy interval, 

5-H ' 



Abstract 



if the produced flux is taken from the Standard Solar Model (SSM). Gal- 
lium detectors, GALLEX, SAGE and GNO, detect additionally pp-neutrinos. 
The pp flux can be calculated subtracting from the gallium signal the rate 
due to boron, beryllium and CNO neutrinos. The ratio of the measured pp- 
neutrino flux to that predicted by the SSM gives the survival probability for 
pp-neutrinos. Comparison with theoretical survival probabilities shows that 
the best (among known models) fit is given by LMA and LOW solutions. 
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The recent measurement of boron v e flux by SNO combined with Super-Kamiokande 
data , gives strong evidence for neutrino oscillations |T[ . It is based on the fact that the v e 
flux measured by SNO through the charged current (CC) interaction u e +d — > e+p+p induces 
in the Super-Kamiokande detector less electrons than observed. The excess of the electrons 
can be produced only by active neutrinos of other flavors: and u T . This flux of active 
neutrinos is thus determined and found to be 3.3a above the zero value. Such analysis was 
suggested earlier in Ref . || , and recently it was further developed in Ref . [Q (for other recent 
analysis of SNO data see [|J1). 

In this Letter we shall demonstrate that the new results obtained by SNO allow us to 
derive the survival probability for boron, beryllium and pp electron neutrinos (for general 
analysis see |7j and for calculation of survival probability for boron neutrinos ||). 

The probability of electron neutrinos to survive on the way from the production point 
inside the Sun to the detection site on the Earth is referred to as survival probability, P ee . In 
case of oscillations, 1 — P ee is the probability of electron-neutrino conversion into neutrinos 
of other flavors. 

The flux of 8 B electron neutrinos determined by SNO via CC-events is $> Ve = 1.75±0.148- 
10 6 cm _2 s _1 |T|, where we used the upper value for systematic error and summed errors 
quadratically. SK measurements provide the flux = & Ue + 0.154($j, m + & Ut ), where 0.154 
is a ratio of cross-sections a v e /a Uee . The comparison of these two fluxes allows us to find 
$ ! , M + $ !/T and $tot = 5V e + 3V M + 5V T , which is equal to 5.44±0. 99 xlO 6 cm~ 2 s _1 Q. Thus the 
survival probability for 8 B electron neutrinos can be found as P ee = $j, e /$tot = 0.32 ±0.065. 
Note that this derivation of P ee does not depend on the SSM flux. The partial oscillation to 
sterile neutrinos, v e — > u s , diminishes further P ee . This possibility is somewhat disfavored 
by the observation that $ tot is close to prediction of the SSM 

The error in the value of P ee indicated above is calculated by adding all errors in quadra- 
ture and it needs further discussion. The uncertainties of and $ tot are correlated 
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and therefore the usual interpretation of the indicated error is allowed only in the limit 
$> Ue <C $tot- I n fact, the 2a- and 3<r-equivalent intervals are asymmetric and greater than 
la error 0.065 multiplied by factor of 2 and 3, respectively. In particular it can be shown 
that 1 — P ee interval has real 3.3cx deflection from zero value. 

The value of P ee is plotted in Figs. 1-3. The calculated value refers to the whole energy 
interval of boron neutrinos measured in Super-Kamiokande: in Figs. 1-3 the horizontal error 
bars show the width of this interval. The value of P ee is plotted in the middle of this interval 
and looks asymmetric in logarithmic scale. As already mentioned above, one should not 
interpret the many standard deviations which separates P ee from unity in the usual way: 
the probability that P ee = 1 corresponds to about 3.3<r and not to ~ 10cr as it might appear 
from the figure. 

The Homestake detector is sensitive to the boron, beryllium and CNO ^-neutrinos. 
Subtracting the contribution of the measured flux of boron i/ e -neutrinos with the standard 
spectrum to the chlorine detector, we determine the contribution of Be+CNO neutrinos to 
the signal. For the production fluxes we use that of the SSM ||: note that the prediction 
for the Be neutrino flux is reliable and the contribution of CNO neutrinos is small. In 
this case the extracted survival probability P ee is not affected by possibility of v e — > v s 
oscillation. The survival probability is plotted in Figs. 1-3 with errors calculated taking into 
account uncertainties in fluxes together with statistical and systematic errors of detection. 
The horizontal error bar refers to the energy interval of Be+CNO neutrinos. The error is 
strongly anticorrelated with the error in the boron flux: a higher (lower) boron flux implies 
a lower (higher) beryllium flux. 

The detected z/ e -neutrino flux of pp-cycle is found from gallium experiments (GALLEX, 
SAGE and GNO), subtracting the contributions from boron, beryllium and CNO neutrinos 
found from the Kamiokande, SNO and gallium experiments. The production flux is taken 
from SSM calculations ||. The calculated pp flux is robust and agrees with the flux found 
independently from solar-luminosity sum rule. The survival probability is plotted in Figs. 1- 
3. The horizontal error bar refers to the pp-neutrino energy interval. The survival probability 
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FIG. 1. Survival probability of electron neutrinos as a function of energy. Data points are 
extracted from the gallium, chlorine and boron-neutrino signals. The horizontal error bars give the 
energy windows of each datum. For the interpretation of vertical error bars see text. The three 
solid curves show the theoretical survival probabilities for the Gribov-Pontecorvo (GP) solution, the 
LOW MSW solution (LOW) and for the large mixing angle MSW (LMA) solution. The solid LMA 
curve corresponds to neutrino production point at the peak of the boron/beryllium production 
zone; the dashed curve - at the peak of the pp production zone. 
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is not affected by oscillation to sterile neutrinos. We shall summarize now the assumptions 
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FIG. 2. Survival probability of electron neutrinos as function of energy. The solid curve shows 
the theoretical survival probability for the small mixing angle (SMA) MSW solution for neutrino 
born at the peak of the boron/beryllium production region; the dashed curve - at the peak of the 
pp production region. 

involved in the calculations of survival probabilities at different energies. For boron neutrinos 
we used only experimental data, neglecting possible oscillation to sterile neutrinos. For the 
other two energy intervals (pp- and Be-neutrinos) we used for total fluxes the production 
fluxes calculated in the SSM, but the calculations of these fluxes are reliable. In addition, we 
assumed that the measured suppression of high energy boron neutrinos (roughly for E > 6 
MeV) is valid also for the lower part of the boron-neutrino spectrum. Since in both cases 
only electron neutrinos produce the signal, the extracted survival probability is not affected 
by oscillation to sterile neutrinos. 

In Fig. 1 the observed survival probabilities are compared with the predictions of the 
LMA MSW, LOW MSW and Gribov-Pontecorvo (GP) [g solutions. For LMA we use one 
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of the best fit solutions from Ref. Q: Am 2 = 3.7 x 10" 5 eV 2 and sin 2 29 = 0.79; for LOW: 



Am 2 = 1.0 x 10~ 7 eV 2 , and sin 2 29 = 0.97 flIQ| , and for the Gribov-Pontecorvo solution 
H sin 2 29 = 1. The LMA and LOW survival probabilities shown are not averaged over the 
production points in the sun. The LOW curve is practically independent of the production 
point. The LMA solid curve is shown for a production point in the boron/beryllium pro- 
duction zone and the dashed curve shows the survival probability for a neutrino produced 
at the peak of the pp region. Averaging over the production region gives a curve close to the 
solid one for boron and beryllium neutrinos, and a curve between the solid and dashed one 
for pp neutrinos (in fact there is little difference between these curves in the energy region 
of pp neutrinos) . 

Only for illustration, we present x 2 /d.o.j '. values for the fits of the data by different oscillation 
solutions: they are 4/3, 5/3 and 10/3 for the LMA, LOW and GP solutions, respectively. 
In fact, this x 2 analysis could be misleading because when the survival probability is used 
as a variable, the probability distribution is not Gaussian. 

In Fig. 2 the observed survival probabilities are compared with those calculated in SMA 
MSW. One of the best fit solutions, with spectral and temporal information included, is given 
by Am 2 = 4.6 x 10~ 6 eV 2 and sin 2 26 = 1.36 x 10~ 3 . The solid curve corresponds again 
to the neutrino produced in the boron/beryllium production region, while the dashed curve 
shows the low-energy part of the survival probability for the neutrino born in the peak of pp 
production region. The confidence level is characterised illustratively by x 2 /d.o.f. « 12/3. 

In Fig. 3 the observed survival probabilities are compared with predictions of the Just- 
So (Am 2 = 4.6 x 10~ 10 eV 2 and sin 2 2d = 0.83) and Just-so 2 (Am 2 = 5.5 x 10~ 12 eV 2 
and sin 2 29 = 0.96) solutions. The qualitatively worse agreement can be characterised by 
X 2 /d.o.f. ~ 11/3 and x 2 /d.o.f. ~ 13/3 for Just-So and Just-So 2 , respectively. 

We repeat again that the quantities x 2 calculated above have only an illustrative char- 
acter: the obtained value of x 2 f° r eac h solution is not connected with probability in a usual 
way. This analysis indicates however the preferred solutions. 

The preferred solutions are LMA and LOW which are characterised by the lowest values 
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FIG. 3. Survival probability of electron neutrinos as function of energy. The solid curve shows 
the theoretical survival probability for the vacuum oscillation solution; for energy below 0.53 MeV, 
only the average probability is shown. The dashed curve is for the Just-So 2 solution. 
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of x 2 - 

In principle the survival probabilities for boron neutrinos can be given for several energy 
bins. Since the observed spectrum is well described by the SSM spectrum, the energy bin 
analysis will further decrease the probability of SMA MSW, which produces a significant 
spectrum distortion in the boron energy window, and will favour the solutions with a flat 
suppression of boron neutrino spectrum, such as the GP, LMA and LOW solutions. Day- 
night dependence will not change the analysis in a significant way, since these solutions have 
little day-night dependence in the boron high-energy window. 

In conclusion, using the data of all solar-neutrino experiments we have derived the 
electron-neutrino survival probabilities. For boron neutrinos only experimental data are 
used, and the SSM flux is not involved. Survival probability decreases in presence of v e — > v s 
oscillation. For Be+CNO and pp neutrinos the SSM calculations are used for production 
(i.e total) fluxes. This is a plausible assumption, because beryllium and pp neutrino fluxes 
are reliably calculated and CNO fluxes are small. Oscillation to sterile neutrinos does not 
affect the extracted survival probabilities for Be+CNO and pp neutrinos. 

LMA and LOW solutions give better fits, in comparison with other solutions. It is 
also possible, in principle, to calculate for boron neutrinos the survival probabilities for 
several energy bins: the likelihood of solutions with an energy independent suppression for 
E > 6 MeV, i.e. LOW, LMA and GP, will increase as compared with the ones that predict 
an energy dependent suppression. 
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